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SUSY-QCD Corrections to Dark Matter Annihilation in the Higgs Funnel 
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We compute the full 0(a s ) SUSY-QCD corrections to dark matter annihilation in the 
Higgs- funnel, resumming potentially large //tan/3 and Ab contributions and keeping all finite 
0(mb, s, 1/ tan 2 (3) terms. We demonstrate numerically that these corrections strongly influence 
the extraction of SUSY mass parameters from cosmological data and must therefore be included in 
common analysis tools such as DarkSUSY or micrOMEGAs. 



INTRODUCTION 

Thanks to the recent mission of the WMAP satellite 
and other cosmological observations, the matter and en- 
ergy decomposition of our Universe is known today with 
unprecedented precision [lj . Direct evidence for the exis- 
tence of Cold Dark Matter (CDM) is accumulating 0,13], 
and its relic density Ocdm can now be constrained to the 
rather narrow range [4| 



0.094 < ft CDM /i < 0.136 



(1) 



at 95% (2(7 ) confidence level. Here, h denotes the present 
Hubble expansion rate Hq in units of 100 km s -1 Mpc -1 . 

Although the nature of Cold Dark Matter still remains 
unknown, it is likely to be composed of Weakly Inter- 
acting Massive Particles (WIMPs), as proposed by var- 
ious extensions of the Standard Model (SM) of particle 
physics. In Supersymmetry (SUSY), a natural candidate 
is the Lightest Supersymmetric Particle (LSP), which is 
stable, if i?-parity is conserved. It is usually the lightest 
of the four neutralinos, denoted Xi or shortly x- 

The Minimal Supersymmetric Standard Model 
(MSSM) depends a priori on 124 soft SUSY-breaking 
parameters, which are often restricted to five universal 
parameters that are imposed at the unification scale 
and can be constrained using data from high-energy 
colliders. As the lightest neutralino relic density depends 
also on these parameters, its computation is another 
powerful tool to put constraints on the parameter space 
and provide complementary information, in particular 
at high energies or masses that would otherwise not be 
accessible at colliders. 

To evaluate the number density n of the relic particle 
with velocity v, one has to solve the Boltzmann equation 



^ = -3Hn - (cr eff v) (n 2 - n\ 



(2) 



with the Hubble rate H and the thermal equilibrium 
density n eq . The present number density uq is directly 
related to the relic density S^cdm^ 2 = ^x n o/Pc °c 
(<7 e ff v}^ 1 , where m x is the LSP mass, p c = 3Hq / (8itGn) 
is the critical density of our Universe, and Gat is the 
gravitational constant || . The effective cross section a c g 



involves all annihilation and co-annihilation processes of 
the relic particle x hrto SM particles, and (cr e ffw) in Eq. 
([2]) signifies the thermal average of its non-relativistic ex- 
pansion (d C c), 

The most important processes contributing to cr e ff are 
those that have two-particle final states and that occur 
at the tree level [6j. Possible final states include those 
with fcrmion-antifcrmion pairs as well as with combina- 
tions of gauge (W ± , Z°) and Higgs (h°, H°, A , 
bosons, depending on the region of the parameter space. 
Processes producing fermions or antifcrmions may be de- 
tectable either directly or through their annihilation into 
photons. In addition, these channels are always open 
(for 6-quarks if m x > 4.5 GeV) in contrast to the other 
channels, which may be suppressed or even closed @, [f| . 

Several public codes perform a calculation of the dark 
matter relic density within supersymmetric models. The 
most developped and most popular ones are DarkSUSY 
[?[ and micrOMEGAs 8]. All relevant processes are imple- 
mented in these codes, but for most of them no (or at 
least not the full) higher order corrections are included. 
However, due to the large magnitude of the strong cou- 
pling constant, QCD and SUSY-QCD corrections are 
bound to affect the annihilation cross section in a signifi- 
cant way. They may even be enhanced logarithmically by 
kinematics or in certain regions of the parameter space. 

In this Letter, we compute these corrections for 
neutralino-pair annihilation into a bottom quark- 
antiquark pair through the s-channel exchange of a pseu- 
doscalar Higgs-boson ^4°. This process dominates in 
the so-called A-funnel region of minimal supergravity 
(mSUGRA) parameter space at large tan j3, which is 
theoretically favored by the unification of Yukawa cou- 
plings in Grand Unified Theories (GUTs) [9j. Supposing 
a WIMP mass of 50-70 GeV, this process has also been 
claimed to be compatible with the gamma-ray excess ob- 
served in all sky directions by the EGRET satellite [loj ]. 
However, the corresponding scenarios may lead to an- 
tiproton overproduction, so that they would not be com- 
patible with the observed antiproton flux 11 1 . 



2 



ANALYTICAL RESULTS 



Denoting the neutralino and 6-quark velocities by 



1 



4m?, 

and /3 b 



1 



4to 2 



(3) 



and the squared total center-of-mass energy by s, the 
properly antisymmetrized neutralino annilation cross sec- 
tion can be written at leading order (LO) of perturbation 
theory as 



1 f3 b N c g 2 Tl u h 



Abb ■ 



2 8ns \ s 



m. 



im A T A \" 



(4) 




a b 




It is proportional to the inverse of the flux factor sv, 
the integrated two-particle phase space s 0b /(8tts), the 
number of quark colors Nq — 3 and the squares of the 
weak coupling constant g, a neutralino mixing factor 

T Alj = ~ (N v - tan OwNij ) (n u cos (3 ~ N 3i sin /?) 

+(*<-► j), (5) 

the bottom-quark mass mi through the Yukawa coupling 
fiAbb — — gm-b tan / (2mw) and the Higgs-boson propa- 
gator. Expanding in powers of it, we obtain ctlo V = 
olo + &lo + O (y i ) with 



flLO = 2&LO — ■ 



N c g 2 Tl n h 



2 

Abb 



47TTO 2 



i ^ — 



'1-5 (6) 



in agreement with Ref. [6(. 

Using standard methods for the virtual one-loop and 
real emission contributions shown on the left-hand side 
of Fig. [1] we compute the 0(a s ) QCD correction in the 
on-shell scheme 
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l + /3b 
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19 + 2/3^ + 3/3^) log 



which agrees with the known result for pseudoscalar 
Higgs-boson decays [ij} and contributes to the total cor- 
rection 



1 



*QCD 



Atop + A 



SUSY 



(8) 



and equivalently for T A - Here, Aq C d = (if-)AqcD 
(^) 2 Ag CD + (^) 3 Ag CD + ... 



FIG. 1: Left: Diagrams for self-energy (top), vertex (center), 
and real (bottom) QCD corrections. Right: Diagrams for 
two-loop top-quark (top) and one-loop self-energy (center) 
and vertex (bottom) SUSY-QCD corrections to the process 
XX -> A -.65. 



In the limit to 2 <C s (0b — > 1), the correction in Eq. 
(|7|) develops a logarithmic mass singularity 



QCD 



a s (s) 



C, 



3, s 9 
"2 l0g ^ + 4 



(9) 



which can be resummed to all orders using the renormal- 
ization group, i.e. by replacing rrib with the running mass 
fhb(s) in the Yukawa coupling h A bb [Hj]- The remaining 
finite QCD-corrections in the MS -scheme are known up 
to 0{a 3 s ) Q 
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*QCD 



a s (s) 



a s (s) 



a s (s) 



(35.94 - 1.36ti/) 



(164.14- 25. 76n/ + 0.259n?). (10) 



A separately gauge-independent 0(a 2 ) correction is in- 
duced by the top-quark loop diagram shown on the upper 
right-hand side of Fig. [TJ Its contribution [15j 
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(11) 

can be important at small values of tan/3, but is largely 
suppressed in the Higgs-funnel region considered here. 

In SUSY, additional 0(a s ) corrections arise through 
the sbottom-gluino exchanges shown in the central and 
lower right-hand side diagrams of Fig. [I] The self-energy 
diagram leads to the mass renormalization Q 



Am b 



a s (s) 



^tan/3)/(m 2 i ,m 2 2 ,m~), 
(12) 



which is for toj, <C tosusy proportional to 
2mi,(Ai > — [i tan 0) 



sin 29h 



to? 

bi 



TO? 
b 2 



(13) 
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i.e. the off-diagonal component of the sbottom mass ma- 
trix, and the 3-point function at zero external momentum 



I(a,b, c) 



ab log f + be log | + ca log | 

(a — 6) (6 — c)(c — a) 



(14) 



In this low-energy (LE) limit, and neglecting Ab with 
respect to the tan /3-enhanced fx, the vertex correction 
equals the mass renormalization up to a factor 
1 / tan 2 j3, so that the total SUSY correction becomes 



A 



(LE) 
SUSY 



a s (s) 



C F 1 
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tan 2 /3 
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(15) 



It has long been known that for large tan/?, Am& can 
be significant and must be resummed by replacing nib — > 
mt/(l -flim^-^ Ami) in the Yukawa coupling /i^bb flit - 
More recently it has been observed that Ab may be of 
similar size as /xtan/3, e.g. in no-mixing scenarios, so 
that its contribution must also be resummed by replacing 
limA b -»o Am 6 -> (linu 6 ^o Am,,)/ (1 + lim M tan/J->o Am 6 ) 
[13]. Our result for the full SUSY-QCD correction A SU sy 
agrees with those in [l8|> [T3| , and we implement the finite 
0{mb 1 s, 1/ tan 2 (3) remainder as described in [171 ] . 



NUMERICAL EVALUATION 

For our numerical study of the impact of QCD, top- 
quark loop, and SUSY-QCD corrections on dark matter 
annihilation in the Higgs-funnel, we place ourselves in a 
minimal supergravity (mSUGRA) scenario with Aq = 
and large tan/3 = 44.5 (54) for fi < (fx > 0), which 
still allows for electroweak symmetry breaking (EWSB) 
in a large region of the scanned m 1 / 2 — mo plane. The 
weak-scale MSSM parameters are then determined with 
SPheno [20| . which includes resummed lim^^o Am,{, cor- 
rections, and the physical Higgs and SUSY masses with 
FeynHiggs [2lj after imposing the current SM masses 
(in particular fhb(mb) — 4.2 GeV and m t = 174.2 GeV), 
gauge couplings (a and sin 2 8w m the improved Born ap- 
proximation, a s (Mz) = 0.1176), and direct and indirect 
SUSY mass limits [22j . For the impact of SUSY spectrum 
calculations on dark matter annihilation see [23| . 

Comparing the observed CDM relic density in Eq. ([T]) 
to the one calculated with DarkSUSY Q, which includes 
the QCD corrections up to C(a 2 ) and where we have 
added the 0(a 3 s ) QCD and 0(a s ) SUSY-QCD correc- 
tions described above, we determine the allowed regions 
in the iriiw — mo plane shown in Fig. [51 The Higgs- 
funnel contribution to a c s rises from 40% for low values 
of mi/2 (or mo for /i < 0) to more than 95%, when m l / 2 
(and m ) is (are) large. For m = 1200 GeV (dashed 
line), the condition (tha — 2iti x )/Ta = is, e.g., satisfied 
when m 1/2 = 840 (970) GeV for /z < (fi > 0), where 
m x ~ 360 (420) GeV. It is obvious from Fig. that the 




600 SOO 1000 

m„ 2 (GeV) 




800 1000 

m„ 2 (GeV) 



FIG. 2: Regions in the mSUGRA mi/2 — mo plane forbidden 
by a charged LSP/no EWSB and favored by the observed 
ficDM h 2 for large tan f3 = 44.5 (54) and fx < (fx > 0). 



LO allowed regions (light) are dramatically changed by 
the 0{a 2 s ) (medium) and 0{a 3 s ) QCD and 0{a s ) SUSY- 
QCD corrections (dark), which reduce <7 e ff by more than 
a factor of two. The increase in Ocdm must therefore be 
compensated by smaller masses. However, Ta is reduced 
by appoximately the same amount, so that on the Higgs 
pole, where Ta is of particular importance, the effect is 
reversed. As expected, the effect is negligible in the focus 
point (very low m^) and co-annihilation (very low mo, 
see also [24|) regions. 

In Fig. [3] we plot the relic density for mo = 1200 GeV 
as a function of m 1( / 2 . Since Ta increases with mx/2, 
in particular for fi > 0, a c g (Oqdm) reaches its maxi- 
mum (minimum) at some distance from the pole (vertical 
dashed line). Further away, the (SUSY-)QCD corrections 
decrease (increase) a e g (^cdm) as expected. Closer to 
the pole, the reduced width becomes important, so that 
the maximum (minimum) approaches the pole. The ef- 
fect of the 0{a 2 s ) QCD corrections already included in 
DarkSUSY and micrOMEGAs [H, [26| is considerably en- 
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FIG. 3: LO (dotted), 0(a 2 s ) (dashed) and 0(a 3 s ) (dot-dashed) 
QCD and 0(a s ) SUSY-QCD (full) corrected predictions for 
ficDM h 2 compared to the experimentally allowed range as a 
function of m 1( / 2 . The position of the Higgs pole is indicated 
by a vertical dashed line. 

hanced by our O(a^) QCD and SUSY-QCD corrections, 
while the top-quark loop contributes less than 0.01% for 
the large values of tan/3 under consideration here. Note 
that the dip at = 420 GeV occurs when m x = mt, 

where the ti annihilation channel is opened |27j | . 

CONCLUSION 

In summary, we have computed the full 0(a s ) SUSY- 
QCD corrections to dark matter annihilation in the 
Higgs- funnel, resumming potentially large /i tan (3 and 
Ai, contributions and keeping all finite 0{mi >1 s, 1/ tan 2 j3) 
terms. We have demonstrated numerically that these cor- 
rections strongly influence the extraction of SUSY mass 
parameters from cosmological data and must therefore 
be included in common analysis tools like DarkSUSY or 



micrOMEGAs. 
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